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This paper consists of two (still only vaguely) related parts: in the first, we briefly review work done 
in the past three years on the "planar equivalence" between a class of non-supersymmetric theories 
_ 1 - (including limiting cases of QCD) and their corresponding supersymmetric "parents"; in the second, 

' we present details of a new formulation of planar quantum mechanics and illustrate its effectiveness 

^ ' in an intriguing supersymmetric example. 

o 

■ 1 Introduction 

o 
\o 
o 



X 



Lattice Gauge Theory (LGT) -an approach where Adriano has distinguished himself 
with many important achievements- is the main tool available at present for extracting 
quantitative results from QCD. This approach has established unambiguously some key 
non-perturbative features of QCD, such as confinement, chiral symmetry breaking, the 
Oh. absence of a [/(l)-Nambu-Goldstone (NG) boson, and more. In spite of these uncon- 

tested successes, LGT remains limited, in its predictive power, by its range of applicabil- 
ity (in terms of number and kind of external particles and of the accessible kinematical 
regions) and by the fact that most of its results, being mere numerical computer outputs, 
often offer limited theoretical insight. 
^ • Few techniques are available today for the analytic study of non-perturbative QCD. 

I Among these, the method of effective chiral lagrangians, based on the general concepts of 

spontaneous symmetry breaking (SSB) and (pseudo) Nambu-Goldstone bosons is one of 
the most widely used. Even though successful, it is also limited in scope to the spectrum 
and interactions of light NG bosons. At the other extreme, heavy quark effective theories 
have been of great value for describing the behaviour of hadrons carrying charm and 
beauty. 

A very different kind of analytic tools is given, finally, by large- expansions of 
various kinds, an approach that is a priori of much wider applicability. Formulated more 
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than three decades ago , large- limits still remain elusive while attracting the attention 

of many theorists^. The relation between the large- A classification of diagrams and loop 
expansions in string theory only enhances the interest in this approach. Unfortunately, 

with the exception of some quantitative results for the rj' mass and interactions?^, large- A^ 
predictions have been mostly qualitative in nature. 

The situation concerning quantitative analytic results in strongly coupled gauge the- 
ories improves considerably if one is willing to consider supersymmetric variants -or 



extensions- of QCD. Predictivity is increased dramatically with the number of super- 
symmetries: while for = 1 this is basically limited to holomorphic quantities, such 
as the superpotential or gauge kinetic terms, for = 2 it extends much further as in 

the famous Seiberg-Witten solutionis. Finally, for = 4, we have the much celebrated 

ADS/CFT correspondence^ between a gauge theory on the boundary and a gravity theory 
in the bulk of anti de-Sitter space. 

In view of the above discussion, one may argue that combining large- ideas and 
supersymmetry should provide as much an analytic predictive power as one may hope 
for. However, both large- and supersymmetry appear to take us further and further away 
from our initial goal: a non-perturbative understanding of (A^ = 3, non-supersymmetric) 
QCD. 

In the first part of this contribution we will summarize work done during the past three 
years on a new large- A^ expansion that has the virtue of connecting (a sector of) ordinary 
QCD to a supersymmetric "cousin", and thus of allowing us to derive both qualitative 
and quantitative properties of the former. In the second part we apply this new large- A^ 
methodology to quantum mechanical gauge systems similar to those obtained by dimen- 
sional reduction of fully space-extended field theories. In particular, we define and solve 
(numerically as well as analytically) a seemingly simple supersymmetric system, which 
exhibits several intriguing features such as a phase transition in the 't Hooft coupling and 
a new form of strong-weak duality. This should provide, eventually, a good and simple 

laboratory for testing the idea of planar equivalence^^. It should also allow an extension 
of recent studies of Supersymmetric Quantum Mechanics!^ to the N ^ oo limit, and to 
make a possible contact with M-theory^. 

In a sense, therefore, the two parts of this contribution to Adriano's fest try to fit with 
the two keywords in its title: Physics and Beauty. Achieving their synthesis represents, 
needless to say, our (yet unattained) final goal. 



2 Large- A^ expansions: Who needs another one? 

As we mentioned in the introduction, very few techniques are currently available for an- 
alytic studies of (non-supersymmetric) gauge theories such as QCD at a non-perturbative 
level. Among the most promising ones, large- A^ expansions play a special role, particu- 
larly because of their conjectured connection to string theories (see e.g.ll^. 

The simplest and oldest large- A^ expansion in QCD is the one suggested in 1974 by 't 

Hooft'^. It considers the limit Nc —>■ (yo while keeping the 't Hooft coupling A = g^N^, as 
well as the number of quark flavours Nf, fixed. Only quenched planar diagrams survive 
to leading order. Quark loops are suppressed by a power of {Nf /Nc) per loop, while non- 
planar diagrams are suppressed by a factor of l/N"^ per handle. This 't Hooft expansion 
led to a number of notable successes in issues such as justifying the validity of the OZI 

rule and in the -q' mass formula^. Unfortunately, nobody succeeded in fully solving QCD 
even to leading order in this expansion. 

For questions where quark loops are important, e.g. in processes with a large number 
of produced hadrons, a better approximation to full QCD is provided by the topological 

expansion (TE)*^, where Nf /Nc, rather than Nf itself, is kept fixed in the large- A^c limit. 
Thus, at leading order, the TE keeps all planar diagrams, including quark loops. This is 



easily seen by slightly modifying 't Hooft's double-line notation by adding a flavour line 
to the single colour line for quarks. In the leading (planar) diagrams the quark loops are 
"empty" inside, since gluons do not couple to flavour. Needless to say, obtaining analytic 
results in the TE is even harder than in the case of the original 't Hooft expansion. 

In Ref.Ha new large- A'^c expansion, which shares some advantages with the TE while 

retaining a significant predictive power, was proposed (for a review, seelH|). Its basic 
idea is actually quite simple. Let us start from ordinary (Nc = 3) QCD with Nf quark 
flavours. Quarks can be described by a Dirac field transforming in the fundamental repre- 
sentation of SU(3)coiour. or, equivalently , in the two-index antisymmetric representation 
(plus their complex conjugates). In extrapolating from A'^c = 3 to arbitrary Nc, the for- 
mer alternative leads to the 't Hooft limit. The new expansion explores instead the latter 
alternative, by representing the quark of a given flavour by a Dirac field in the two-index 
antisymmetric representation. In that case, taking the Nc ^ oo limit at g^Nc and Nf 
fixed does not decouple the quark loops since, for large Nc, the number of degrees of 
freedom in the antisymmetric field scales like the one in the adjoint, i.e. as N^. For rea- 

sons explained in"^, it has been referred to as the orientifold large- A'^c limit. The leading 
order of this new expansion corresponds to the sum of all planar diagrams, in the same 
way as in TE, but with the crucial difference that quark loops are now "filled", because 
the second line in the fermion propagator is also, now, a colour line. 

The orientifold large- A'^c limit is, therefore, unquenched. Its 't Hooft-notation dia- 
grams look, modulo reversal of some arrows in the fermion loops, precisely as those of 
an SU(A'^c) gauge theory with A'^^ Majorana fields in the adjoint representation, a theory 
that we may call adjoint QCD. Adjoint QCD can be seen, in turn, as a softly broken ver- 
sion of supersymmetric Yang-Mills theory (SYM) with Nf — 1 additional adjoint chiral 
superfields. The soft-breaking terms are just large mass terms for the scalar fields that 
make them decouple. 



2. 1 Proofs of planar equivalence 

The planar equivalence between M = 1 SYM theory and orientifold field theories 
amounts to the following statement: an SU(A'^c) gauge theory with two Weyl fermions 
in the two-index antisymmetric representation (i.e. one Dirac antisymmetric fermion) is 
equivalent, in a bosonic subsector and as A'^c ^ oo , to = 1 gluodynamics. 

In Ref.l2], a perturbative proof of the planar equivalence was provided and a non- 

rrri 

perturbative extension was outlined. In a subsequent paper a refined non-perturbative 
proof of planar equivalence (extended to A'^^ > 1) was given. Its basic idea is the com- 
parison of generating functionals of appropriate gauge-invariant correlators in the parent 
and daughter theories. This is done by, first integrating out their respective fermions in 
a fixed gauge background and, then, averaging over the gauge field itself. The first step 
produces a fermionic determinant, which, in turn, can be expanded as a sum of Wilson 
loops computed for different representations. Using some group-theory identities as well 
as certain factorization properties of Wilson loops at large A^, it was possible to relate 
the generating functionals of the parent and daughter theories for a subset of suitably 

identified sources. Recently, a lattice version of the non-perturbative proof of has been 

given , together with suggestions on how to check the equivalence by realistic lattice 
simulations. 

It is worth noting that the two theories are not fully identical. In particular the colour- 



singlet spectrum of the orientifold theory consists only of bosons and does not include 
composite fermions at Nc ^ oo. 



2.2 SUSY relics in QCD 

We will now use planar equivalence to make predictions!^ for one-flavour QCD, keeping 
in mind that they are expected to be valid up to corrections of the order of = 1/3 
(barring large numerical coefficients): 

(i) Confinement with a mass gap. Here we assume that large- A'^c J\f = 1 gluodynamics 
is a confining theory with a mass gap. Alternatively, if we start from the statement that 
one-flavour QCD confines, we arrive at the statement that M = I SYM theory shares 
that property, while the mass gaps are dynamically generated in both theories. 

(ii) Degeneracy in the colour-singlet bosonic spectrum. Even/odd parity mesons (typ- 
ically mixtures of fermionic and gluonic colour-singlet states) are expected to be degen- 
erate. In particular, 
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mz., 

^ = 1 + 0{1/Nc) , one-flavour QCD , (1) 

where rj' and a stand for 0^ and 0+ mesons, respectively. This follows from the exact 
degeneracy in = 1 SYM theory. Note that the a meson is stable in this theory, as there 
are no pions. This prediction should be taken with care (i.e. a rather large numerical 
coefficient in front of l/Nc may occur), since the t]' mass is given by the anomaly (the 

WV formula'^), whereas the a mass is more "dynamical." The degeneracy between even 
and odd-parity mesons should improve at higher levels on the expected Regge trajectory. 
In order to check this relation in real QCD, let us assume that the a mass is not very 

sensitive to the number of flavours. On the other hand, according to the WV formula™ 
and neglecting quark masses, the r]' mass scales like y/Nj; we can therefore extrapolate 
relation ([T) to obtain a prediction for real QCD 

m^' ~ VSnia. (2) 

Although the a bump is very broad, it is encouraging that the above relation is indeed in 
qualitative agreement with the position of the enhancement in the appropriate vrvr channel. 

(iii) Bifermion condensate. M = I SU(A'^r;) gluodynamics has abifermion condensate 

[El that can take Nc distinct values: 

(AA)fe ~ M^^ ^-r/N,^2i^k/N, ^ ^^3 ^i{9+2nk)/N, ^ k = 0, I, . . . , - 1, (3) 

with 

r = — id 

9^ 

and c a calculable numerical coefficient. The finite- A'c orientifold field theory is non- 
supersymmetric, and here we expect (taking account of pre-asymptotic 1/Nc corrections) 
Nc — 2 degenerate vacua with 

{^L^R)k' ~ exp <^ — + i- 



<72(Afe + 4/9) iV, 



c 



~ dt^ ^^P jV - 2 J ' ^ =0,l,...,iVc-3. (4) 

The term 4/9 in (IHl is due to the one-loop /? function of the orientifold field theory, 
h = 3Nc + |, while A'^c — 2 in (IDl is twice the dual Coxeter number of the antisymmetric 
representation (fixing the coefficient of the axial anomaly). Finally, c' is a normaUzation 
factor that we will discuss below. 



2.3 Analytic estimate of the quark condensate in QCD 

The relation between the quark condensate in one-flavour QCD and the gluino condensate 

in SYM theory can be pushed further by appealing to the fact that the quark condensate 
of QCDoR must agree with the gluino condensate at = oo and must vanish at = 2. 
Thus, its value at = 3 (which is nothing but the quark condensate of one-flavour QCD) 
can be obtained by interpolating between these two values. The final outcome of such an 

analysis is the following analytic formula: 
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= -5^m'(A0.)) " '■■ exp^-^Jk(l/3) . (5) 

where /j, is the renormalization scale for the condensate, (3q and (3i are the one- and two- 
loop /3-function coefficients, 7 is the one-loop anomalous dimension of the condensate, 
and k{l/N) stems for further, unidentified 1/A^ corrections. 

In a recent papeJ^an attempt was made to extend these considerations to an arbitrar 



number of quark flavours. The idea (actually first proposed for different reasons in^S 
is to add to QCDqr rif quarks in the (traditional) fundamental representation of SU(A'^), 
so that they become irrelevant in the large- A^ limit. This new class of theories has been 
dubbed QCD^r' . At A^ = 3, such a theory is nothing but QCD with Nj flavours where 
A'^^ = Uf + I. One has to determine which sectors of QCDqr' can be mapped into 
corresponding quantities of SYM theory, in particular since the former theory has a NG 

sector which is absent in SYM. When this is done^^^one obtains an interesting general- 
ization of dS]) to Nf > 1 and can even argue that the connections should be particularly 
accurate for the "realistic" case of A'^^ = 3. For three-flavour QCD the relevant values 
are po = 9, /3i = 32, 7 = 4. 

In order to calculate the condensate ([5]l we need to know the value of the 't Hooft 
coupling A at a scale /i that we choose to be ^ = 2 GeV. The Particle Data GroupES 
quotes as (2 GeV) = 0.31 ± 0.01, which corresponds to A(2 GeV) = 0.148 ± 0.010. 
We therefore choose to plot the function 

(^^')2GeV 3 3 44 / 1\ 

— ;2^X 27 exp — -- , (6) 



in a range of A. Comparison between the above theoretical prediction and present deter- 
minations of the condensate is shown in Fig. 1 . Clearly the result supports the validity of 
planar equivalence within the expected precision at A^ = 3. 
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Figure 1: The quark condensate expressed as — (y MeV)^ as a function of the 't Hooft coupling A. The zhlcr 
range of the coupling, 0.138 < A < 0.158, and the lattice estimate -(259 ± 27 MeV)^ define the shaded 

region. 

3 Planar Quantum Mechanics 

In this second part of the paper, we turn from QFT to the simpler case of Quantum 
Mechanics and try to define a "Planar Quantum Mechanics" that hopefully mimicks the 
large- limit of a system whose degrees of freedom are x (bosonic or fermionic) 
matrices. The hope is that, eventually, this technique can be extended to QFT and provide 
a new handle on solving interesting gauge theories -such as YM, SYM and QCD- in the 
large- limit. It may also provide useful information on whether (and how) "planar 
equivalence" may or may not work. 

3. 1 Fock states, and matrix elements at large N 

Our quantum mechanical system is defined by a Hamiltonian H, which is polynomial in 
A'^ X A^-matrix annihilation and in creation operators aij, a^ij. Symmetry under the U {N) 
transformations is ensured by taking H to be the trace of such a polynomial. Creation 
and annihilation operators satisfy standard commutation relations 



[aij,a'ki] — ^ii^jk 



(7) 



We shall work in the eigenbasis of the U(A'')-invariant number operator B = Tr{a^a). 
Therefore basis states are constructed by acting on the Fock vacuum with the invari- 
ant building blocks, or "bricks", Tr[(a^)"], n = 1,2, . . . , N. Following the Cayley- 
Hamilton theorem there are N independent bricks for the U(A^) group. 

A general state is a product of all possible bricks and their powers. It will prove useful 
to define a cut Hilbert space Hb, i.e. a space restricted to states with no more than B 
bosonic quanta. Such a space is spanned by all polynomials of bricks with maximum 
B-th order. Obviously there are many states with a fixed number of bosonic quanta n. 

In the large- limit, things simplifyl^ni considerably: for given n, there is only one 
relevant state. It is created by the single-trace operator: 

\n) = ^Triia^rm , (8) 

where Mn is a suitable normalization factor. All states that are created by products of 
traces are non-leading in the sense that they give rise to non-leading matrix elements. A 
second simplification is that leading operators also have a single trace form. These two 
basic rules of the Planar Calculus are best illustrated by the explicit examples below. 
Normalization of planar states 
The normalization factor AA„ in ^ reads 

= (0|rr[a"]rr[(at)"]|0) = (0|(12)(23)...(nl)[l'2'][2'3'] . . . [n'l']|0), (9) 

where we have simplified the notation, e.g. ai^^i^ (12) and a^i7,ig [7'8'] and 
all indices are summed. Using d?]) it is easy to see that the maximal power of is 
achieved by contracting indices between () and [] starting from the middle and working 
outwards. This produces the famous planar diagrams. One way to picture it is to represent 
()()...() as a circle with n pairs of short lines, representing indices, pointing inside and 
[][]...[] as a smaller circle with lines pointing outside. Leading contributions result when 
opposite pairs of lines match: this gives n loops, i.e. a factor N'"-. Because of the cyclic 
properties of a trace in one gets n such contributions. Moreover, according to Q, 
each contraction gives a factor 1. The final result reads 

l^ivw., „>0, 

I 1, n = u. 

Example of matrix elements 

As another example, we calculate the matrix element of a typical term in a generic Hamil- 
tonian: 

Hn+2,n = 9^{n + 2\Tr[a)a)a)a]\n). (11) 

Consider a state 

Tr [a^a)a)a] Tr [(a^ )"] 1 0) = [12] [23] [34] (41) [l' 2'] [2'3'] . . . [n' l'] 1 0) 

= n[l'2][23][32'][2'3'][3'4']...[n'l']|0), (12) 



n equal terms result from commuting one annihilator, (41), all the way until it hits the 
vacuum state. Each contraction (|7]) gives just a factor 1. Pictorially, one can attach a 
small circle with 3 creators and one annihilator to a big circle with n pairs of indices 



in n equivalent ways. The resulting state is proportional to |n + 2). Collecting the 
normalization factors gives 



y/{n + 2)/ nN 



i^n+2,n = g^n-V(n + 2|Tr[(Qt)("+2)]|0)=g2n - x 



(n + 2|rr[(a^)("+2)]|0) = g'^N^n[n + 2), n > 0. (13) 

Indeed, such an interaction term depends only on A = g^N, showing the relevance of the 
t'Hooft coupling. To complete this example, we also calculate a matrix element of the 
conjugate operator 

Hn-2,n = 9^{n - 2\Tr[a^aaa]\n) = -^(n - 2|Tr[a1'aaa]Tr[(a^)"]|0). (14) 

The leading contribution comes from the three subsequent contractions of adjacent a' s 
and s, which give a factor A^^. Again this can be done in n ways. Taking into account 
the renormalization Mn J^{n-2) gives finally 



n—2,n 



n > 2, (15) 



which again shows the t'Hooft scaling and, upon the replacement n —>■ n + 2, confirms 
the hermiticity of the planar hamiltonian. 



3.2 An anharmonic oscillator 

Let us now apply planar rules to find the spectrum of the C/(oo) -invariant anharmonic 
oscillator. The Hamiltonian 

H = Ti[p^]+g^TT[x'^]=T + V (16) 

is bounded from below and represents a perfectly well defined system per se. On the 
other hand, there exists a whole family of similar Hamiltonians (supersymmetric or not), 
which have been obtained from gauge-invariant field theories by reducing the system to a 
single point in space. For example, reduced two-dimensional Yang-Mills gluodynamics 
is described solely by the kinetic term of Eq. (fT6l ) with a global U{N) invariance as the 
reminder of the local gauge symmetry of the original, space-extended theory. 
As a first step, introduce matrix creation and annihilation operators 

x = ^{a + a)), p=-^Xa-a}), (17) 

and rewrite each term of H in the normal-ordered form: 

2 2 ^ 
+ (^-]^ + Ng'^^{a^'^ + a^) + {l + 2g'^N)o)a 

+ ^(atVa^)+5VV + <7'(at'a + ata3) (18) 



A note on the normal ordering may be useful here. For higher powers of a's and a^'s the 
normal ordering may spoil the group structure. For example 

: aiko) kjajia\i := kjO^ uaiuaju (19) 
is not a trace, hence it does not contribute to the planar limit. On the other hand 

: a'ikakjajio'ii := a) ika^iakjaju (20) 

can be brought to the trace form - it preserves the group structure - and consequently 
contributes to the leading behavior. As an example consider one of the quartic terms in 
the expansion of the potential in Eq.([T6l). 

(2,2) = Trya^a^"^ + aa) aa) + aa) a) a + aaa) + aa^ a + a^"^ a^\. 



Commuting all annihilation operators to the right and retaining only single traces gives 

(2,2) = 2N'^ + mTr[a)a]+ATr[a)'^a\ (21) 



Out of the six quartic terms only four preserve the trace structure after bringing them to 
the normal ordered form. The remaining two do not, hence they are non-leading. 

The calculation of all matrix elements of ([18]) closely follows the earlier examples. 
After some algebra one obtains 



< n\H\n > 




(l + 3A)n, n>2, 
l + 2A)n, n<2, 



, -i + 2A ^n{n + 2), n > 0, 
<n + 2\H\n> =<n\H\n + 2> \ (22) 



2 



i + A ^/2N 



n 



. A\TT\ \TT\ . A I TAA/n(n + 4), n > 0, 

<n + A\H\n> =<n\H\n + 4> = <^ 4 v v 

I 2-^^ n = 0. 

This result illustrates an important feature of all large- calculations (also valid in QFT). 
Namely, not all matrix elements scale with A. There exists a class of "superleading" 
contributions which are divergent in the 't Hooft limit. They result from vacuum diagrams 
and should be treated separately!. At the moment, we just neglect them, i.e. we use in 
the following the subtracted Hamiltonian H — N'^{1 + A)/2 and ignore the first row and 
column where the non-scaling vacuum diagrams contribute only. 

We may now go ahead and diagonalize the so modified Hamiltonian H. To this end 
we introduce a cut-off n < B, find numerically the spectrum of H, and increase B until 

the convergence is reached!^. 

Our results are shown in Fig. |2]and can be summarized as follows. 

• Satisfactory convergence of lower eigenenergies, with increasing cut-off, is achieved 
for A < ~ 4. In this region the system resembles an effective harmonic oscillator 
with almost (but not exactly) equidistant "infinite volume" levels. This seems to 

be a general feature of the large- dynamics!^ 



"We are grateful to E. Onofri for useful discussions on this issue. 
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Figure 2: Cut-off dependence of the spectrum of the anharmonic oscillator in the even-parity sector. 



• With A approaching Ac ~ 4 the mass gap vanishes and all levels fall towards 
zero with increasing cut-off. Such a behaviour has been studied before and is the 

characteristic sign of the continuous spectrum in the infinite-volume limil)231. 

• For A > Ac all levels fall with the cut-off towards —oo;H appears to be unbounded 
from below. 

The last point may seem unnatural, signalling some trouble for the planar approach. In 
fact it provides one more argument that a quantum mechanics turns into a field theory 
in the large- limit. Namely, bound states with negative energies are of course quite 
common in both systems. However only a field theory, with its multiparticle sectors, can 
account for an infinite series of negative states^ 

Nevertheless this emergence of the negative eigenenergies provides an important les- 
son (and a warning) about possible intricacies of the planar limit. The original Hamilto- 
nian (fT6l ) is perfectly well defined and positive for any finite N. This raises the question 
of the relation between the "true" large- limit, defined as solving non-perturbatively a 
system at finite N and then taking N large, and the planar limit discussed here. There 
may be many reasons for the non-commutativity of the two, the most obvious one being 
the ad hoc subtraction of superleading contributions. The model discussed hereafter does 
not have this deficiency. 



3.3 A very symmetric supersymmetric system 

We now add fermionic degrees of freedom described by matrix creation and annihilation 

operators f\ f, {fijfli} = SuSjk ■ To avoid vacuum diagrams, we require that super- 
symmtery charges (hence also a Hamiltonian) annihilate the empty Fock state. A simple 



choice 

Q = Tr[fa\l+ga% = Tr[f\l + ga)a] H = {Q,Q^} = Hb + Hp, (23) 
or exphcitly (a trace is always implied) 

J. . J. 2 4- o, O 4-2 o 

Hb = a^a + g{a^ a + a'a )+g a' a , 

Hp = Pf + g{Pfia^ + a) + f\a^ + a)f) (24) 

+ g^f^afa^ + f^aa^f + f^fa^a + f^a)fa). 

This Hamiltonian conserves the fermion number F = f. In each fermionic sector the 
planar basis {\F,n)} is now created by the single trace with F fermionic and n bosonic 
creation operators. It is now a simple, but somewhat tedious, exercise to calculate the 
matrix elements of H in this representation. We obtain for the first two fermionic sectors: 

{0,n\H\0,n) = (1 + A)n - A(5„i, 



{0,n + l\H\0,n) = {0,n\H\0,n + l) = Vx^n{n + 1), (25) 

{l,n\H\l,n) = (1 + A)(n + 1) + A, 

{l,n + l\H\l,n) = {l,n\H\l,n + l) = VX{2 + n). (26) 

The numerical calculation of the spectrum proceeds in the same way as for the anhar- 
monic oscillator (see Fig.©. 



'JW thanks P. Menotti for a very instructive discussion on this subject. 
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Figure 3: Cut-off dependence of the spectrum of H, in the F = sector, in a range of values for A. 
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Figure 4: First 10 energy levels of _ff , in _F = and F = 1 sectors, at A = 0.5. 



Now, however, the salient features of this system are: 

• The spectrum is non-negative for all values of the 't Hooft coupling. This confirms 
the suspicion of vacuum diagrams being responsible for the earlier problems. 

• For A away from 1 the eigenenergies (at infinite cut-off) are again almost equally 
separated with a common mass gap depending on A. 

• There is a phase transition at Ac = 1. At this point, and only there, the spectrum 
loses its mass gap and becomes continuous. 

• In the vicinity of Ac we observe a critical slowing down: higher cut-offs are re- 
quired to achieve satisfactory convergence. 

• There is an excellent boson-fermion degeneracy, see Fig. ID Supersymmetry is 
unbroken at infinite cut-off but is quite badly broken at finite cut-off near the critical 
point Ac. 

• In both phases there is an unpaired bosonic SUSY vacuum with zero energy. This 
is the empty Fock state, which, by construction, is annihilated by our Hamiltonian. 

• While moving across the transition point, at finite cut-off, members of supermul- 
tiplets rearrange; see Fig. \5\ In particular a second ground state with zero energy 
appears in the strong coupling phase. 

• Some symmetry between the strong and weak coupling regions emerges. 

It turns out that the problem can be solved analytically, offering better understanding 
of these features. We shall now discuss each one of them in turn. 



Supersymmetry 

With the rules of Sect.2 the matrix representation of supersymmetry charges is easily 
obtained in the planar limit: 

l|Qt|o,n) = \/n, {l,n - 2\Q^0,n) = 0,n) = b^/^. (27) 

Obviously the generator has matrix elements only between F = and F = I sectors. 
One can now readily check the sum rule 

^(0,n'|Q|l,m)(l,m|Q^|0,n) = {0,n'\H\0,n) , (28) 

m 

which follows from H = {Q,Q^. In our representation its test becomes a simple 
exercise in matrix multiplication. Notice that even the exceptional diagonal contribution 
at n = 1 is correctly reproduced by the anticommutator of SUSY charges. 
A Uttle more involved test employs the identity 

(0,n'|//2|0,n) = {0,n'\QHQ^\0,n), (29) 

which again turns into straightforward matrix algebra. Notice that the sum rule resulting 
from ( [29l) relates matrix elements of H in different fermionic sectors. 

We conclude that the planar approximation preserves supersymmetry. On the other 
hand the cut-off in terms of the number of bosonic quanta explicitly breaks SUSY. Conse- 
quently the fermion-boson degeneracy is observed numerically only at sufficiently large 

cut-offs. It is conceivable that SUSY invariant cut-offs can be found^S. 



The second ground state 

To understand the existence of the second massless state, let us introduce the composite 
creation and annihilation operators and a„, which create the planar basis in the new 
Hilbert space: 

|0,n) = a^„|0). (30) 
In terms of these operators the F = Hamiltonian reads 

oo / oo \ 

ff{F=o) ^ ^t^^^ + ^ „(i + 62)at„a„ + 6^n(n + l)a'^„a„+i + h.c. . (31) 

n=l \n=l / 

Then one can formally construct the following state 

I0h = i:(li) ^10).. (32) 

which is annihilated by Q^, Eq. (|27] ). hence also by the Hamiltonian (1251 ). It is normaliz- 
able only for A > 1, thereby explaining the puzzle found numerically. The e merg ence of 

this second ground state in the strong coupling phase causes the Witten indexl^to jump 
by one unit (in the sectors with F = 0, 1) across the phase transition. According to the 
Feynman-Kac relation, the thermal partition function also reveals such a discontinuity 
at zero temperature. However, since higher levels collapse to zero at Ac, there is also a 
(5-like contribution at this point: c.f. Fig.[6l 
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Figure 5: Rearrangement of the F = and F = 1 spectra, around A = 1, for increasing cut-offs. 




Figure 6: A dependence of the Witten index and the partition function, at /3 = 6, for different cut-offs. 



Strong-weak duality 

The qualitative symmetry of the energy levels between the strong and weak coupling 
phases, seen in Fig. [51 finds its explanation in a rather intriguing exact duality between 
the two regimes. Consider the Hamiltonian ( |26l ) in the F = 1 sector. It is evident that the 
"reduced" Hamiltonian 

H = ^{H-b^) (33) 
is symmetric under the replacement b ^ 1/6. Therefore the eigenenergies of H satisfy 

The same reciprocity relation in the F = sector does not hold at the operator level, 
c.f. (1251) . On the other hand, because of the supersymmetry, the duality works there as 
well. However, since there exists a second bosonic vacuum, in this sector duality relates 
eigenenergies with different indices: 

i(£;f=o)(6)-62)=6(<T°^(l/6)-l), b<l. (34) 
The analytic solution 

Finally, and rather surprisingly, it turns out that the Hamiltonian (1251 ) can be analytically 

diagonalized. More details can be found irJ^, here we outline only the main steps of the 
derivation. 

1. The F = Hamiltonian is strikingly simple when written in terms of the linear 
combinations of a„: 

oo 

H = ^ B^Bn , Bn = Vnan + bVn + lun+i , 

n=l 



2. The states created by the B}^ operators 

\Bn)=Bl\0) = ^/^\n) + bV^^\n + 1) . (35) 

form a non-orthonormal basis. We shall find the eigenstates IV') of the reduced Hamil- 
tonian by expanding them into the basis and constructing the generating function 
f{x) for the expansion coefficients: 

oo oo 

IV') = Cn\Bn) ^ fix) = J2 Cn^"- (36) 
n=0 n=0 

3. The simple way in which H acts on the \Bn) basis implies the following first-order 
differential equation for f{x) 

H\i;)=e\i>) ^ w{x)f'{x)+xf{x)-bf{0)-f'{0) = ef{x), (37) 

where w{x) = (x + b){x + 1/b), and the eigenvalues of H are simply related to those of 
H: E = 6(e + b). 

4. A solution of (|37] ) is straightforward (see againl^for details). The final expression 
for the generating function reads 

where F{a, b, c; x) is the standard geometric function, and the quantization condition, 

(e + 6)/(0) = , (39) 

determines the discrete series of eigenenergies in the F = sector. 

This analytic solution nicely explains all results discussed above. In particular: 

• The absolute value of /(O) reveals a series of zeros in a, c.f. Fig. |7l which agree 
with the large cut-off limit of the numerical eivenvalues of ( |25] ). 



• The duality among massive eigenvalues follows immediately from ( |39b and the 
symmetry of ( [38] ) under the substitution b l/h and a — > 1 — a. 

• The collapse of the eigenenergies to zero at the critical point is evident from E = 
a{b'^ — 1), and so is the fact that the roots are bounded by nearby poles of the 
/3 function. 

• The generating function of the second ground state can be easily obtained by set- 
ting a = in Eq. ( [38] ) for 6 > 1. This gives 



which indeed corresponds to the expansion ([32]) . On the other hand, one cannot set 
a = in the 6 < 1 solution - the second vacuum does not exist for 6 < 1. 



b = 2 . 




Figure 7: The quantization condition as a function of a. The first four zeros are clearly visible. To see higher 
zeros the a resolution of the plot must be increased. 



Given (l38l ) and (|39l ) one can study the flow of the eigenvalues in both phases. Perhaps 
the most interesting is the behaviour in the vicinity of the critical point. From the known 
asymptotics of the hypergeometric functions around A = 1 we conclude that the first 
eigenvalue tends to zero as —(1 — A)/log(l — A) when A 1^, which results in a 
vanishing first -and infinite second- derivative. 

One can also quantitatively study the free energy near the phase transition, which 

appears to be stronger than in the Gross-Witten modeM 



4 Higher-F sectors 



The structure of higher-F sectors, rather than being a simple repetition of the F = 0, 1 
sectors, appears to exhibit novel interesting features. This already follows quite directly 
from looking at the free theory (A = 0). Even in that case the structure of the multiplets, 
and the way they should rearrange into supermultiplets, is highly non-trivial. This is best 
illustrated in a plot that resembles the old Chew-Frautschi plot (CFP), in which angular 
momentum is plotted against squared mas^ 

As one can see, while at all mass levels there is just one F = and one F = \ state 
that naturally form SUSY doublets, this is not the case for the F = 2, F = 3 states. 
Their number, at high enough occupation numbers, is different, the excess always being 
in favour of the F = 3 sector. However, supersymmetry must be satisfied, and this is 



Linearity of the (Regge) trajectories in this plot was perhaps the first evidence for a string-like structure 
of hadrons, with the inverse of the Regge slope reinterpreted as the string tension. 




:. 1 1 1 1 III 



3 4 
E 



Figure 8: Chew-Frautschi plot of the F = and F = 1 eigenenergies at weak 't Hooft coupling. 



achieved non-trivially by having some of the F = 3 states pair with those with F = A. 
Once more, this leaves some of the F = 4 states unpaired but, continuing the exercise 
upwards in the CFP, shows that, eventually, every state finds its supersymmetric partner. 

An amusing mathematical digression can be made at this point: our states are, in many 
ways, what mathematicians would call "binary necklaces" because they have a cyclic 
structure and are made out of just two kinds of "beads", a bosonic and a fermionic one. 
One immediately finds on the simplest examples that, in general, necklaces of a given 
(even) length do not split in an equal number of necklaces with even or odd fermions 
(there are always more bosonic than fermionic necklaces). This would immediately vio- 
late SUSY, except that Pauli's principle comes to the rescue by killing a subset of bosonic 
necklaces. Such a result, which can obviously be formulated as a purely mathematical 
combinatorial problem, appears to be newB 

Coming back to our CFP, we may discuss further which linear combinations of the 
degenerate states of a given F and E form SUSY doublets. The task is facilitated by 

introducing another conserved quantityl^^ 



C^iQlQ], [C,H]=0 



(41) 



and by noting that = H^. Clearly the eigenstates can be classified according to their 
"C-parity" i.e. according to whether C = ±E. We may even consider the combination 
of C and fermion number Cp = (— 1)^C. States with C = +1 (—1) are clearly those 
annihilated by (Q). Most of the states turn out to have Cp = — 1 (all of them in 
the F = 0, 1 sectors), but there also are some "unnatural" ones with Cp = +1, which 
are very important for the full matching of bosons and fermions discussed above. While 
the analysis of generic values of F needs a more systematic -and probably computer- 
based- approach, we have been able to deal by brute force with the F = 2, 3 sectors 

confirming all the results discussed above. There is, however, another surprise: As 



■^GV wishes to thank Professor J.C. Yoccoz for this private information. 



A ^ 1, all states again go to ii^ = 0; at A > 1, two states with F = 2 stay at £^ = 0, 
while all other states move up to positive E. The question of whether new E = states 
keep popping up at large A in even higher-F sectors is currently being explored. 

5 Outlook 

We would like to conclude this paper by outlining a possible programme for connecting 
its two parts in a reasonably near future. 

A first step would consist in generalizing our SQM model to the case in which bosons 
and fermions belong to different representations of the symmetry group U{N). In the 
case of bosons in the adjoint representation and fermions in the two-index (anti)symmetric 
representation (plus its complex conjugate), it would be interesting to check whether (at 
least parts of) the bosonic spectra of the latter theories coincide with those of the super- 
symmetric case in the planar limit. 

Next, we would like to extend our planar calculus to systems with rotational sym- 
metry, such as those obtained by dimensional reduction of (supersymmetric) Yang-Mills 
theories in = 3 + 1 and D = 9 + 1 dimensions. This basically requires generalization 
of the planar rules for more than one species of bosons and fermions. A good step in 
this direction is being made by studying our supersymmetric quantum mechanical model 

with an arbitrary number of fermions^. 

Should these first two steps succeed, the next one would be to extend the whole ap- 
proach to the planar Hamiltonians of full fledged quantum field theories, with or without 
supersymmetry, and in particular to the orientifold large- limit of QCD described in the 
first part of this paper. 

Finally, comparing the N = oo results with those that can be obtained at A^ = 3 by 

the method of Ref. ^ ^ , would give an estimate of how good the orientifold 1/A^ expansion 
is at A^ = 3 and of the extent to which supersymmetric gauge theories can make accurate 
predictions for actual QCD. 
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